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We calibrate the efFective-one-body (EOB) model to an accurate numerical simulation of an equal- 
mass, non-spinning binary black-hole coalescence produced by the Caltech-Cornell collaboration. 
Ahgning the EOB and numerical waveforms at low frequency over a time interval of ~ lOOOM, and 
taking into account the uncertainties in the numerical simulation, we investigate the significance and 
degeneracy of the EOB adjustable parameters during inspiral, plunge and merger, and determine 
the minimum number of EOB adjustable parameters that achieves phase and amplitude agreements 
on the order of the numerical error. We find that phase and fractional amplitude differences be- 
tween the numerical and EOB values of the dominant gravitational wave mode /i22 can be reduced 
to 0.02 radians and 2%, respectively, until a time 20M before merger, and to 0.04 radians and 7%, 
respectively, at a time 20M after merger (during ringdown). Using LIGO, Enhanced LIGO and 
Advanced LIGO noise curves, we find that the overlap between the EOB and the numerical ^22, 
maximized only over the initial phase and time of arrival, is larger than 0.999 for equal-mass binary 
black holes with total mass SO-ISOM©. In addition to the leading gravitational mode (2, 2), we com- 
pare the dominant subleading modes (4, 4) and (3, 2) for the inspiral and find phase and amplitude 
differences on the order of the numerical error. We also determine the mass-ratio dependence of 
one of the EOB adjustable parameters by calibrating to numerical inspiral waveforms for black-hole 
binaries with mass ratios 2:1 and 3:1. The results presented in this paper improve and extend recent 
successful attempts aimed at providing gravitational-wave data analysts the best analytical EOB 
model capable of interpolating accurate numerical simulations. 

PACS numbers: 04.25.D-, 04.25. dg, 04.25. Nx, 04.30.-w 



I. INTRODUCTION 

The first-generation gravitational-wave detectors — 
the Laser Interferometer Gravitational Wave Observa- 
tory (LIGO) [1, 2], GEO [3] and Virgo [4] — have op- 
erated at design sensitivity for a few years, providing 
new upper limits for several astrophysical sources. They 
are now undergoing an upgrade to Enhanced LIGO and 
Virgo"*"; this will improve their sensitivity by a factor of 

2. The second-generation interferometers, Advanced 
LIGO [5] and Advanced Virgo, will start operating in 
2013-2015 with an overall improvement in sensitivity by 
a factor of ~ 10, thus increasing the event rates for many 
astrophysical sources by a factor of one thousand. 

One of the most promising sources for these detectors 
is the inspiral and merger of compact binary systems of 
black holes. The search for gravitational waves (GWs) 
from coalescing binaries and the extraction of parameters 
are based on the matched- filtering technique [6, 7], which 
requires a rather accurate knowledge of the waveform 
of the incoming signal [8]. In particular, the detection 
and subsequent data analysis of GW signals are made by 
using a bank of templates modeling the GWs emitted by 
the source. 

The efFective-one-body (EOB) formalism was intro- 
duced [9, 10] as a promising approach to describe an- 
alytically the inspiral, merger, and ringdown waveforms 
emitted during a binary merger. Necessary inputs for 



the EOB approach include high-order post-Newtonian 
(PN) results [11] for two-body conservative dynamics, 
radiation-reaction force, and gravitational waveforms. 
For compact bodies, the PN approximation is essentially 
an expansion in the characteristic orbital velocity v/c or, 
cquivalently, in the gravitational potential, GM / {rc^), 
with r the typical separation and M the total binary 
mass. The EOB approach, however, docs not use the 
PN results in their original Taylor- expanded forms (i.e., 
as polynomials in v/c), but instead in some resummed 
forms [9, 10, 12-17]. The latter are designed to incor- 
porate some of the expected non-perturbative features of 
the exact results. 

As it is now possible to produce very accurate numeri- 
cal simulations of comparable mass binary black-hole co- 
alescences (see e.g. [18-25]), we can compare in detail the 
EOB predictions with numerical results, and when nec- 
essary, introduce new features into the EOB model in or- 
der to improve its agreement with the numerical results. 
This is an important avenue to LIGO, GEO and Virgo 
template construction, as eventually thousands of wave- 
form templates may be needed to detect the GW signal 
within the detector noise, and to extract astrophysical in- 
formation from the observed waveform. Given the high 
computational cost of running the numerical simulations, 
template construction is currently an impossible demand 
for numerical relativity alone. 

This paper builds upon a rather successful recent ef- 
fort [22, 26-31] aimed at producing the best analytical 



EOB model able to interpolate accurate numerical sim- 
ulations. Other approaches based on phenomenological 
waveforms have also been proposed [32, 33]. Here we 
calibrate the EOB model to the most accurate numer- 
ical simulation to date of an equal-mass, non-spinning 
binary black-hole merger, that has been produced with 
a pseudospectral code by the Caltcch-Cornell collabora- 
tion [21, 23]. Taking into account the uncertainties in 
the numerical simulation, we investigate the significance 
and degeneracy of the EOB adjustable parameters and 
determine the minimal number of adjustable parameters 
that achieves as good agreement as possible between the 
numerical and EOB GW's phase and amplitude. In ad- 
dition to the leading GW mode {i,m) = (2,2), we also 
compare the leading subdominant modes (4, 4) and (3, 2). 
By reducing the phase difference between the EOB and 
numerical inspiral waveforms of black-hole binaries with 
mass ratios q = mi:m.2 of 2:1 and 3:1, we explore the 
dependence of one of the adjustable parameters on the 
symmetric mass ratio v =■ mi ra-ij {mi + 7712)^. 

The paper is organized as follows. In Sec. II, we briefly 
review the EOB dynamics and waveforms. In Sec. Ill, 
we calibrate the EOB model to the numerical simulation 
of an equal-mass non-spinning binary black-hole coales- 
cence and determine the region of the parameter space 
of the EOB adjustable parameters that leads to the best 
agreement with the numerical results. We also discuss 
the impact of our results on data analysis, and calibrate 
the EOB model with inspiral waveforms from accurate 
numerical simulations of non-spinning black hole bina- 
ries with mass ratios 2:1 and 3:1. Sec. IV summarizes 
our main conclusions. Finally, the Appendix compares 
the numerical him extracted with the Regge-Wheeler- 
Zerilli (RWZ) formalism with the him obtained by two 
time integrals of the Newman-Penrose (NP) scalar ^['4'". 



II. EFFECTIVE-ONE-BODY MODEL 

In this section we briefly review the EOB dynamics 
and waveforms, focusing mainly on the adjustable pa- 
rameters. More details can be found in Refs. [9, 10, 13, 
15, 22, 26, 28-31]. Here we follow Refs. [22, 28]. 



toise radial coordinate r. 
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In terms of Pr, the non-spinning EOB Hamiltonian is [9] 
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with the effective Hamiltonian [9, 13, 34] 
H''''{r,pr,,p^)=^iH''' 



Pr+Mr) 
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The Taylor-approximants to the coefficients A{r) and 
D{r) can be written as [9, 13] 
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The functions A{r), D{r), Ak{r) and Dk{r) all depend on 
the symmetric mass ratio v through the i/-dependent co- 
efficients ai{v) and di{v). These coefficients are currently 
known through 3PN order (i.e. up to fc = 4) and can be 
read from Eqs. (47) and (48) in Ref. [22]. Previous inves- 
tigations [15, 22, 28-31] have demonstrated that, during 
the last stages of inspiral and plunge, the EOB dynamics 
can be adjusted closer to the numerical simulations by in- 
cluding in the radial potential A(r) a pseudo 4PN (p4PN) 
coefficient a^{v). This coefficient has so far been treated 
as a linear function in v, i.e. a^{v) = XqI', with Ao a 
constant^. In this paper, however, we shall also explore 
the possibility of going beyond this linear dependence, 
such that 



A. Effective-one-body dynamics 

We set M = mi -I- 7712, fj. = mi m2/M = v M , and use 
natural imits G = c = I. In absence of spins, the motion 
is constrained to a plane. Introducing polar coordinates 
(r, $) and their conjugate momenta {pr,p^), the EOB 
effective metric takes the form [9] 

dsls = -A{r) df + ^dr^+r^ (dO^ +sm^ 6 d^^^ . (1) 



05 (i^) = i^(Ao + Ai , (7) 

where Aq and Ai are constants. In order to assure the 
presence of a horizon in the effective metric (1), a zero 
needs to be factored out from A{r). This is obtained by 
applying a Pade resummation [13]. The Pade coefficients 
for the expansion of A{r) and D{r) at p4PN order are 
denoted A\{r) and Dl{r), and their explicit form can be 
read from Eqs. (54) and (59) in Ref. [22]. 



Following Ref. [16, 34], we replace the radial momentum 
Pr with pr_^ , the conjugate momentum to the EOB tor- 



1 Note that Ao was denoted A in Ref. [28] , and as in Refs. [22, 29- 
31], 
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The EOB Hamilton equations are written in terms of 
the reduced (i.e., dimensionlcss) quantities H'^°^^ [defined 
in Eq. (3)], f = t/M, and h^flM [10]: 
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with the definition O = d^/dt. Furthermore, for the <!> 
component of the radiation-reaction force we use the non- 
Keplerian Padc-approximant to the energy flux [12, 35] 
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where wn = O^^'^, V$ = flrQ, and = r [^{;{r,p^)Y^^ . 
Here ip is defined by Eqs. (66)-(68) of Ref. [22]. As the 
EOB Hamiltonian is a deformation of the Schwarzschild 
Hamiltonian, the exact Keplerian relation = 1 

holds. The quantity in Eq. (12) is given by Eqs. (39) 
and (40) in Ref. [22]^ and it uses the Taylor-expanded 
energy flux (as given by Eq. (19) in Ref. [22]) in the form 
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where we combine the known test-mass-limit terms [36] 
with a p4PN adjustable parameter As [22].^ 

The radial component of the radiation-reaction force 
Trir,Pr,,P<i>) in Eq. (10) was neglected in previous stud- 
ies [22, 26, 28-31] because Ref. [10] showed that for 
quasi-circular motion, in some gauges, it can be set 
to zero. Furthermore, it was shown in Ref. [10] that 
if the motion remains quasi-circular even during the 
plunge, J-'r{r,pr,,p<i,) does not affect the dynamics con- 
siderably. However, since we are trying to capture effects 
in the numerical simulations which go beyond the quasi- 
circular motion assumption, we find it interesting to add 



^ Note that here we use the Pade approximants with factorized 
logarithms, as originally proposed in Ref. [12], but wo set dlso = 
1, so that the GW energy flux depends only on the two adjustable 
parameters Vp^ic a-^d Ag. 

3 Note that in Ref. [22] the p4PN contribution in the GW energy 
flux also included the term uBg log vq. Since we found apprecia- 
ble degeneracy between As and Sg, we disregard Bs, i.e., we set 
Bs = 0. 



Tr{r,Pr,,P't>) [see Eq. (3.18) of Ref. [10] and the discus- 
sion around it] . We set 



^r{r, Pr.-.P'i.) 



where a^p {v) is an adjustable parameter. 

Finally, the tangential force described by Eq. (12) ap- 
plies only to quasi-circular motion. This tangential force 
could also in principle contain terms describing the de- 
parture from quasi-circular motion during the last stages 
of inspiral and plunge. There are several ways to include 
such non-quasi-circular (NQC) terms [15, 22, 31]; here 
we do so by replacing the quantity on the right-hand 
side of Eq. (11) [but not the T,^ on the right-hand side 
of Eq. (14)] with NQ"^.?*, where 



(15) 



and {v) is an additional adjustable parameter. The 
form of this NQC correction will be discussed further 
in Sec. IIIB. Note that alternative NQC terms have 
been proposed in the literature — for example, in Ref. [15] 
the authors used p^/(p$/r^) while Ref. [31] employed 
pI /(rfi)^. In summary, in the notation of Ref. [22], the 
EOB model used here is ^'^F^/H^ with adjustable pa- 
rameters {05 (z^) , -Wpoio (i^) , a^g (z/) , a^^ {v),Ag,}. 



B. EOB waveform: Inspiral & Plunge 

Having the inspiral dynamics in hand, we need to com- 
pute the gravitational waveform him- Reference [21] 
compared the numerically extracted gravitational wave- 
form /122 to the FN result with amplitude expressed as 
a Tai/Zor-expansion [37, 38]; even when expanded to 
3PN order, the amplitude disagreed by about one per- 
cent at times several hundred M before merger. As pre- 
vious investigations [29-31] have shown, more accurate 
agreement with the numerical /122 amplitude can be ob- 
tained by applying several resummations to the Taylor- 
expanded /i22 amplitude. These resummations have re- 
cently been improved using results in the quasi-circular 
test-particle limit [17]. We follow Ref. [17] and write the 
EOB modes /i£,„ as 
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where R is the luminosity distance from the binary, and 
with 
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where ffeff and Joff are effective sources that in the test- 
particle, circular-motion limit contain a pole at the EOB 
light ring (photon orbit); here HeS is given in Eq. (4), and 
JcS = Pis>vn is equal to the orbital angular momentum p$ 
normalized to the circular-orbit Newtonian angular mo- 
mentum v^^. The quantities Tgm, Sg„i, Pim, P'tm '^^^^ 
read from Eqs. (19), (20), (23), (25), (CI), (C4) and (C6) 
in Ref. [17], respectively. More specifically, T^™ is a re- 
summed version [16] of an infinite number of leading log- 
arithms entering the tail effects; Sem is a supplementary 
phase [16] which corrects the phase effects not included in 
the complex tail factor; p^m and pf^ are the resummed 
expressions of higher-order PN effects as recently pro- 
posed in Ref. [17] in the test-particle circular-orbit limit. 
The latter resummation was proposed to cure, among 
other effects, the linear growth with £ of the IPN correc- 
tions in the Taylor-expanded amplitude. 

Furthermore, motivated by the PN expansion for 
generic orbits, to include NQC effects in him we write 
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As we shall discuss in detail below, for the (2,2) mode, 
one of the four adjustable parameters a^^= in Eq. (18) 
will be fixed by requiring that the peak of the EOB /122 
occurs at the same time as the peak of the EOB orbital 
frequency [31] (i.e., at the EOB light-ring); this requires 
no matching to a numerical waveform. Another of the 
a^^^ will be fixed by requiring that the peak amplitude of 
the EOB and numerical waveforms agree. The final three 
a^^^ parameters will be determined by minimizing the 
overall amplitude difference with respect to the numerical 
waveform. We note that an alternative NQC factor has 
been proposed in Ref. [31], notably 1 -\- a p^J {Vi^ 
We shall compare those different choices below. 



C. EOB waveform: Merger & Ringdown 

The merger-ringdown waveform in the EOB approach 
is built as follows [10, 26, 28, 30, 31, 35]. For each mode 
(£, to) we write 



N-l 



merger— RD 



(19) 



where n is the overtone number of the Kerr quasi-normal 
mode (QNM), N is the number of overtones included in 
our model, and Atmn are complex amplitudes to be de- 
termined by a matching procedure described below. The 
quantity aimn = ^imn — ioiimm where the oscillation fre- 
quencies LOimn > and the inverse decay-times aimn > 0, 



are numbers associated with each QNM. The complex 
frequencies are known functions of the final black-hole 
mass and spin and can be found in Ref. [39]. The final 
black-hole masses and spins can be obtained from sev- 
eral fitting formulae to numerical results [28, 34, 40, 41]. 
Here we use the more accurate final black-hole mass and 
spin computed in Ref. [23]: Mbh/M = 0.95162±0.00002, 
a/il/BH = 0.68646 ±0.00004. While these numbers differ 
from the predictions of the fitting formulae in Ref. [28] 
by only 0.3%, such disagreement would be noticeable in 
our comparison. The matching time imatch(^) 
justable parameter that will be chosen to be very close 
to the EOB light-ring [10] when matching the mode ft.22. 

The complex amplitudes Aimn in Eq. (19) are deter- 
mined by matching the EOB merger-ringdown waveform 
with the EOB inspiral-plunge waveform. In order to do 
this, N independent complex equations are needed. In 
Refs. [10, 26, 28, 35, 42], the N equations were obtained 
at the matching time by imposing continuity of the wave- 
form and its time derivatives 

^ I insp— plunge/ ,£m \ 7 merger— RD /,fm \ 

^-^k im \^ match )— ^^k irn I'match/ ' 

(fc = 0,l,2,..- ,7V-1), (20) 

and we denote this approach point matching. In Refs. [30, 
31], the comb matching approach was introduced. In 
this approach, N equations are obtained at N points 
evenly sampled in a small time interval At^^^^^ centered 
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Finally, the full (inspiral-plunge-mcrgcr-ringdown) EOB 
waveform reads 



h _ I insp-plunge /)/,£m 
IHrn — ll-im "If-match^ 



-t) + h 
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(22) 

The point matching approach gives better smoothness 
around the matching time, but it is not very stable nu- 
merically when iV is large and higher order numerical 
derivatives are needed. As we include eight QNMs in our 
ringdown waveforms, we find that the comb matching 
approach is more stable. To improve the smoothness of 
the comb matching we use here a hybrid comb matching: 
We choose a time interval At^^^j^ ending at i^tch' 
impose the continuity of the waveform at — 4 points 
evenly sampled from t^^tch " ^Catch to Citeh' but we 
also require continuity of the first and second order time 



derivatives of the waveform at tf^t^h" 



A-iim „ J j-im 
^'•match thatch' 



thus guaranteeing the continuity of him- Furthermore, 
we fix t^tch to be the time when the EOB orbital fre- 
quency reaches its maximum, and tune Ai^^^^^j^ in the 
range 2.5M-3.5M depending on the EOB dynamics. 
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It is worth noticing that the lowest frequency among 
the eight QNMs included in our merger-ringdown wave- 
form is MllI227 ~ 0.44, which is larger than the EOB 
inspiral-plunge waveform frequency MLLi{t'^^^^^J ~ 0.36. 
Therefore, generically the EOB GW frequency will grow 
very rapidly from Mcu ^ 0.36 to Mu! ^ 0.44 immediately 
after the matching time, and this growth can be much 
more rapid than what is seen in the numerical simula- 
tion. We find that we can avoid this rapid growth by 
carefully fine-tuning the matching interval Ai^^^^^j^, and 
this is what we do for the comparisons presented here. 
Quite interestingly, we find that the matching can 
be made much less sensitive to At'^^^^^^ if we include a 
pseudo QNM that has a frequency Muj{t'^^^^yJ 0.36 
and a decay time comparable to that of the highest over- 
tone T227 ~ 0.7M. We refer to this QNM as pseudo be- 
cause its frequency and decay time do not coincide with 
any of the QNMs of our final Kerr BH [39, 43]. Although 
we do not use this pseudo QNM in the present analy- 
sis, we expect that its inclusion can help when matching 
higher modes of equal and unequal mass binaries and we 
shall consider it in the future. 



III. CALIBRATING THE 
EFFECTIVE-ONE-BODY WAVEFORMS TO 
NUMERICAL RELATIVITY SIMULATIONS 

We shall now calibrate the EOB model against a nu- 
merical simulation of an equal-mass non-spinning bi- 
nary black hole. This simulation was presented as run 
"30cl/N6" in Scheel et al. [23], and the inspiral part 
of the waveform was used in previous comparisons with 
PN models [21, 22]. In addition to the NP scalars ^-l™ 
extracted from this simulation, we will be using grav- 
itational waveforms him extracted with the RWZ for- 
malism [44-47]. The Appendix discusses details of the 
numerical implementation used to obtain him from the 
RWZ scalars, and presents a comprehensive comparison 
of the numerical and RWZ /i£,„ waveforms. Consis- 
tency between the two wave-extraction schemes is good, 
with phase differences less than 0.02 radians for the (2,2)~ 
mode until about a time 20M after the peak of |ft.22|- 

Because we have more experience with the NP scalars 
during the inspiral, and because ^'4^ appears to behave 
better than RWZ h22 during ringdown (see Fig. 14 in 
the Appendix), we prefer to use the numerical data. 
Therefore, during the inspiral phase, we will calibrate 
the EOB adjustable parameters by comparing the sec- 
ond time derivative of EOB /122 against the numerical 
^'4^. During the plunge- merger phase, when the time 
derivatives of the waveform vary most rapidly, it is more 
difficult to calibrate the EOB /i22 since the resummation 
techniques in the EOB model were aimed at providing us 
with the best /i22. Therefore, around time of merger, we 
shall calibrate the EOB h22 to the RWZ /i22- Note also 
that data analysis is based on him, further motivating our 
choice to build the best EOB model for him- Neverthe- 



less, after calibration, in Sec. IIIC, we show comparisons 
of the EOB waveforms with both the numerical RWZ /i22 
and ^'1^. The ringdown part of the numerical waveform 
is not used in the calibration of the EOB parameters; the 
QNMs are determined solely from the mass and spin of 
the final hole. 



A. Waveform alignment and uncertainties in 
numerical waveforms 

As previous investigations [21, 26, 31, 48, 49] have 
shown, the phase error between two waveforms depends 
crucially on the procedure used to align them in time 
and phase. For the inspiral phase, we shall adopt here 
the alignment procedure introduced in Ref. [22] (see also 
Rcf. [33]) that consists of minimizing the quantity 

S(At, A<j>)= r [<t>i{t) -Mt- Ai) - A(t>f dt , (23) 
Jti 

over a time shift At and a phase shift A<f>, where (f>i{t) 
and (/)2(i) are the phases of the two waveforms. This 
alignment procedure has the advantage of averaging over 
the numerical noise and residual eccentricity when align- 
ing numerical and EOB waveforms. The range of inte- 
gration (ii, is chosen to be as early as possible, where 
we expect the PN-bascd EOB waveform to be most valid, 
but late enough so that it is not contaminated by the junk 
radiation present in the numerical initial data. Moreover, 
the range of integration should be large enough for the 
integral to average over noise and residual eccentricity. 
Here we fix ti = 1040M and t2 = 2260M (measured from 
the start of the numerical waveform), so that we include 
three full cycles of phase oscillations due to eccentricity. 

Using this alignment procedure, we estimate the er- 
rors on the numerical by comparing ^f^^ computed 
at different numerical resolutions and/or using different 
extrapolation procedures. In particular. Fig. 1 summa- 
rizes the phase errors for a set of numerical ^'4^ computed 
in Ref. [23]. The numerical waveform labeled "N6, n=3" 
(identical to the run "30cl/N6, n=3" from [23]) is the 
reference numerical waveform used throughout this pa- 
per unless otherwise noted. This waveform is the most 
accurate waveform from Ref. [23], extracted at various 
radii and then extrapolated to infinity. The waveforms 
with different values of n vary the order of the extrap- 
olation and are used to quantify the uncertainty in the 
phase due to extrapolation, while those labeled by N5 
(as opposed to N6) are from a simulation with a lower 
numerical resolution and are used to quantify the uncer- 
tainty due to numerical truncation errors. Figure 1 also 
includes a comparison between waveforms extracted at 
finite coordinate radius Tex = 225M. 

Extrapolation with n = 2 leads to systematic errors 
in the extrapolated waveform (see, Fig. 10 of Ref. [21]), 
which in turn results in a systematic error in At. There- 
fore, the blue dashed line in Fig. 1 represents a possibly 
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FIG. 1: Numerical error estimates. Phase difference 
between numerical waveforms, wiien aligned using the 
same procedure as employed for the EOB-NR alignment [see 
Eq. (23)]. "N6" and "N5" denote the highest- and next-to- 
highest numerical resolution, n denotes the order of extrapo- 
lation to infinite extraction radius, and "r = 225 M" denotes 
waves extracted at finite radius r = 225 M. The data are 
smoothed with a rectangular window of width lOM; the light 
grey dots represent the unsmoothed data for the N5-N6 com- 
parison at Tex = 225Af . 



overly conservative error estimate. The feature of the 
solid brown curve around t ~ 3700M is due to an issue 
with data processing of the lower resolution 'N5' run. 

The primary use of Fig. 1 is to assess numerical errors 
relevant for the calibration of the EOB inspiral phase. 
By construction of the alignment procedure, this figure 
shows the numerical errors for waveforms that are aligned 
in the interval [ii,i2], several orbits before merger. Cal- 
ibrating the EOB inspiral phase in this manner is ap- 
propriate, because it ensures that early in the inspiral, 
the EOB-modcl and the numerical simulation agree well, 
i.e. that wc expect little de-phasing at lower frequencies. 
This is important for waveform templates of low mass bi- 
naries, where the early inspiral waveform lies in LIGO's 
sensitive frequency band. 

Figure 1 shows that the numerical ^'4^ waveforms are 
accurate to a few hundredths of a radian until very close 
to merger, when compared with our alignment proce- 
dure. Furthermore, Fig. 15 in the Appendix demon- 
strates that NP and RWZ waveforms differ by only 0.02 
radians through inspiral and merger. Therefore, wc shall 
adopt a deviation of 0.02 radians between EOB- and NR 
inspiral-wavcforms as our goal for the EOB inspiral cali- 
bration. The horizontal line in Fig. 1 indicates this phase 
difference of 0.02 radians and it will be our requirement 
when calibrating the EOB values of ^4^. The numeri- 
cal phase errors exceed 0.02 radians at times t = 3660A/, 
3850A/, 3900A/ and 3933M, respectively, and so our goal 
will be for EOB to agree to 0.02 radians at least up to 



EOB-dynamics 
adjustable parameters 


EOB-waveform 
adjustable parameters 


t'pole(i') 

aRRi") or a^l,{u) 
Aa 


a^"Mz = l,...4 



TABLE I: Summary of all possible adjustable parameters of 
the EOB model considered in this paper. As we shall discuss 
in the main text, we will not need all of these parameters. In 
particular, we find that for the black-hole binary simulations 
investigated here, the choices a^J^(j/) = = a^g^(z^),Ag = 0, 
^match(i^) s-t peak of the EOB orbital frequency, allow 
the numerical and EOB values of the GW phase and ampli- 
tude to agree within numerical error. Furthermore, we find 
that for an equal-mass black-hole binary coalescence it is suf- 
ficient to set 05(1/) — vXo [see Eq. (7) with Ai = ] and 
calibrate Ao, Upoie(l/4), Af?^2,t,i,(l/4) and a,''==(l/4). For an 
equal-mass black-hole binary coalescence it is even possible 
to calibrate only one EOB-dynamics adjustable parameter, 
Ao [see Eq. (7)] and let Wpoie — > 00. Finally, for an unequal- 
mass binary inspiral it is sufficient either to set Ai = 0, use 
the value of Ao from the equal-mass binary case, and calibrate 
Vpoia{i^)', or alternatively to let Wpoie —> 00 and calibrate both 
Ao and Ai in 05(1/) [see Eq. (7)]. 



t « 3900 Af. The choice of 0.02 radians is motivated by 
the goal of bringing the disagreement between the EOB 
and numerical phases at least to the level of the numerical 
error (sec Fig. 15). 



B. Tuning the adjustable parameters of the 
equal-mass efFective-one-body dynamics 

Wc divide the adjustable parameters into two groups 
and tune them separately in two steps. The 
first group of FiOB- dynamics parameters includes 
{a5{iy),Vpoioiiy),a^%^{v),a^^{iy), As}. These parameters 
determine the inspiral and plunge dynamics of the EOB 
model and affect the mergcr-ringdown waveform only 
indirectly through the waveform's phase and frequency 
around the matching point. [Wc note that the inspi- 
ral phase is independent of the parameters Qi'^, see 
Eq. (18).] These parameters are calibrated to the nu- 
merical NP \l/4^. The second group of FiOB-waveform 
parameters includes {af^-, i^^^i^, and Ai^^^i^}, and af- 
fect only the plunge-merger-ringdown but not the inspi- 
ral EOB waveform. These parameters are calibrated to 
the numerical RWZ ft,22- AH the possible adjustable pa- 
rameters of the EOB model employed in this paper are 
summarized in Table I. In the first step of our calibration 
procedure, we reduce the phase difference before merger 
by tuning the EOB-dynamics parameters. In the sec- 
ond step, we use these fixed values of the EOB-dynamics 
parameters, and tune the EOB-waveform parameters. 

Among the EOB-dynamics parameters, 05 (^) and 
'ypoie(i') are the most important as they affect the en- 
tire quasi-circular evolution of the inspiral. The two 



7 



radiation-reaction parameters a^^^ and are intro- 
duced to adjust the dynamics of late inspiral when we 
expect that the quasi-circular assumption is no longer 
valid. The p4PN parameters in the energy flux, Ag, also 
influences the entire evolution, but we find that Ag, is 
strongly degenerate with 05(1/4) throughout the inspi- 
ral until a time ~ lOOM before merger. Based on these 
considerations, we shall tunc 05(1/4) and Upoic first and 
consider o^^(l/4), 0^^(1/4) and Ag only when exploring 
how to further improve the late evolution. 

Therefore, in our first step, we set o^p^(l/4) = 
aRR(V4) ^ As ^ and vary 05(1/4) and ?;poic(f/4). 
Applying the alignment procedure presented at the be- 
ginning of Sec. Ill A, we shift each EOB '^^^ in time and 
phase to agree with the reference numerical waveform at 
low frequency, and determine the time when the phase 
difference between the numerical and EOB wave- 
forms becomes larger than 0.02 radians. We denote this 
reference time as trei- 

Figure 2 is a contour plot of the time t^^f in the 
a5(l/4)-Wpoie(l/4) parameter space. For all points in- 
side the largest contours (blue curves), the associated 
EOB 5*1^ phase evolutions agree with the numerical ones 
up to t = 3660Af, which is the earliest reference time 
considered in Sec. Ill A. In order to get EOB models 
that have phase differences less than 0.02 radians until 
t = 3900 A/, 05(1/4) and i'poie(l/4) have to be inside 
the innermost two separate thin contours (red curves). 
One might view these contours as encompassing all val- 
ues of 05(1/4) and Wpoic(l/4) that are consistent with 
the numerical inspiral waveform, given the fixed choices 
of the various other EOB parameters. There are 05(1/4) 
and Wpoio(l/4) values that make the EOB phase differ- 
ences less than 0.03 radians until t = 3933M, but not 
less than 0.02 radians until t = 3933M (the latest refer- 
ence time). We find that phase errors of the EOB 5*4^ 
corresponding to the upper left contours in Fig. 2 grow 
rapidly after t = 3900Af, whereas phase errors of EOB 
corresponding to the lower right contours grow only 
mildly until around t — 3940Af . For this reason, we shall 
restrict the tuning of the other adjustable parameters 
to the lower right region of Fig. 2 inside the innermost 
contour. As a reference set, we choose 05(1/4) = 6.344 
and i'poie(l/4) = O.SS."' We note that the latter value 
is rather different from the value obtained in Ref. [31] 
when 05(1/4) = 6.25 is used. This is due to differences 



* We note that in Ref. [28], the authors suggested as best value 
05(1/4) = 15. However, the EOB model used in Ref. [28] differs 
from the one employed in this paper, the main difference being 
the GW energy flux. More importantly, the procedure used in 
Ref. [28] to calibrate 05(1/4) was different. It was based on max- 
imized overlaps with white noise. The best value for 05(1/4) was 
obtained by requiring large overlaps, say > 0.0975, for several 
mass ratios and {£,m) modes (see Fig. 2 in Ref. [28]). Finally, 
the accuracy of the numerical waveforms employed in this paper 
differ from the ones in Ref. [28]. 
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FIG. 2: In the parameter space of the EOB-dynamics ad- 
justable parameters as (1/4) and Upo1g(1/4) we show the con- 
tours of the time fref at which the phase difference between 
the numerical "30cl/N6, n=3" and EOB ^'4^ becomes larger 
than 0.02 radians. Note that the innermost red contours cover 
two disjoint regions. The inset shows the effect of numerical 
uncertainty: The filled contours are the t^-^f = 3850Af and 
3900Af contours from the main panel. The open contours are 
identical, except computed using the "30cl/N6, n=2" numer- 
ical ^f^^. The reference model is shown as a black dot. 

between the EOB models — for example Ref. [31] em- 
ploys the Pade-resummed GW energy flux with constant 
logarithms, whereas we use the Pade-resummed GW en- 
ergy flux with factorized logarithms. 

Quite interestingly, looking more closely at the red 
lines in the right corner of Fig. 2, as Vpoie increases, we 
find another possible reference set 05(1/4) = 4.19 and 
''^poic — *■ 00. With this choice, the pole in the Pade flux 
of Eq. (12) disappears. 

In order to understand whether further tunings of 
radiation-reaction effects by adjusting the parameters 
(wpoic, a^p, Opp, Ag) can modify the phasing during 
plunge, we compute how sensitive the phasing is to 
radiation-reaction effects once the binary has passed the 
last stable orbit (LSO) defined as {dH"^ /dr)i,so = = 
(d'^H^^ /dr'^)i^so. Reference [50] pointed out that the 
phasing during the plunge is not affected much by radia- 
tion reaction, but driven mostly by the conservative dy- 
namics. We want to quantify the latter statement more 
fully. 

In order to do this, wc need to define when the plunge 
starts. In the absence of radiation reaction, the plunge 
starts beyond the LSO where r = rj^so, ^ = ^hso and 
P$ = P^^'^ ■ But in the presence of radiation reaction, 
Ref. [10] observed that there is not a unique <lso at which 
the conditions r = tlsOj — ^^lso and p$ = p]^'~' are 
satisfied. In fact, the above conditions may happen at 
different times (see Fig. 12 in Ref. [10]). Indeed, for 
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FIG. 3: For the case a5(l/4) = 6.344 and Upoio(l/4) = 0.85 
(^8 ~ 0, — and a^j^ = 0), we show the phase difference 
between the numerical and FOB mode /i22 versus the numer- 
ical GW frequency Muj22 for FOB models in which the GW 
energy flux is shut down at several FOB orbital frequencies. 
The vertical line marks the maximum FOB orbital frequency. 



the case 05(1/4) = 6.344 and Wpoio(l/4) = 0.85, we find 
that with radiation reaction, r{t^gQ) ~ ri^so, "^(^lso) ~ 
c^LSO and p^itllo) = p^so ^^crc t^^^ = 3914.50M, 
i^gQ = 3919.83M and i[go = 3885. 53M, and where the 
orbital frequencies, corresponding to the three different 
^Lso values are Mfl = 0.975,0.106, and 0.074, respec- 
tively. Following Ref. [10], we will say that the plunge 
starts during the time interval spanned by the values of 
^LSO which in this case is ^lso ~ 34M before merger. 

In Fig. 3, we show the phase difference between the 
numerical and FOB /122 as a function of the numerical 
GW frequency Ma'22 for FOB models in which the GW 
energy flux is suddenly shut down at several FOB or- 
bital frequencies. The cyan curve in Fig. 3 is obtained 
when the GW energy flux is not shut down. Note that 
in this case the phase difference increases fast close to 
the FOB matching point, which is marked by the ver- 
tical line in Fig. 3. The phase difference can change 
considerably its shape (including the sign of the slope 
close to the FOB matching point) when the energy flux 
is shut down before Mfl = 0.12-0.13, but it docs not 
change much, especially the fast increase close to the 
matching point, when the energy flux is shut down af- 
ter Mil ~ 0.12-0.13, immediately after the LSO defined 
by the condition '.^(^lso) ^ "^lso above. 

This study suggests that it is difficult to modify the 
behaviour of the FOB phasing during plunge by tuning 
only the adjustable parameters entering the radiation- 



reaction terms or the GW energy flux, a;^J^(z/), a^^{iy) 
and Ag, i'poie(i^)- The behaviour of the FOB phasing 
during plunge is more sensitive to adjustable parameters 
in the FOB conservative dynamics, e.g., 05(1/) at 4PN 




0.8 

Vpole(l/4) 



FIG. 4: Fffect of Oj^j^ on contours of acceptable FOB param- 
eters. The solid contours are the trcf = 3850M and 3900M 
contours from Fig. 2. The open contours shifted to the lower- 
right 

are the same, but computed with — 0.5 instead of 
'^RR ~ 0. The reference model is shown as a black dot. 



order or ae{i') at 5PN order, etc. However, the param- 
eters ai{i/) also affect the phasing during the very long 
inspiral, and a careful tuning is needed to reach excellent 
agreement both during inspiral and plunge. 

Nevertheless, it is possible to modify the behaviour 
of the FOB phasing during the late inspiral by tuning 
As, a^p^(l/4) and a5^^(l/4) together with 05(1/4) and 
'ypoio(l/4). As an example of this, we redo the contour 
plot shown in Fig. 2, but with a^f^(l/4) = 0.5 instead 
of zero. The result is shown as dashed curves in Fig. 4. 
We still find FOB models that have phase differences 
less than 0.02 radians until t = 3900Af . In particular, 
with the reference value Wpoic(l/4) = 0.85 and choosing 
05(1/4) = 6.013, we find that the behaviour of the FOB 
phasing is substantially modified only for the last 40M 
of evolution before merger. In this case, the change in 
phase difference is in the range of 0.01-0.1 radians, and 
the slope of phase difference at the matching point can 
change sign. Similar results are obtained when repeating 
this analysis with app (l/4) or As different from zero. We 
also observe that the effect on the dynamics of the ad- 
justable parameter 0^^(1/4) is almost equivalent to the 
effect of the adjustable parameter app(l/4), except for a 
minus sign and a different scaling. So it is not necessary 
to consider both of these radiation-reaction adjustable 
parameters. 

Although time consuming, in principle it is possible 
to perform a comprehensive search over the complete 
set of the FOB-dynamics parameters 05(1^), fpoie(i^), ^81 
a^^iv) or (v). However, at this point there is no need 
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FIG. 5: We compare the numerical and EOB h^i amplitudes 
when the EOB model with reference values as (1/4) = 6.344 
and iipoic(l/4) = 0.85 are used. We show the EOB ampli- 
tudes without the NQC corrections and the EOB amplitude 
with the NQC terms suggested in Ref. [31], where the NQC 
parameters take the values a — 0.75 and e = 0.09. When the 
NQC corrections are not included, we show the EOB ampli- 
tude of Eq. (16a) which uses the resummation procedure of 
Ref. [17], and also the EOB amplitudes of Eq. (16a) when the 
Pade-resummations Pi and P3 of P22 suggested in Ref. [17] 
are applied. Note that in this plot, the EOB amplitudes do 
not contain the merger-ringdown contribution. 



to further improve the EOB evolution close to merger, 
and achieve better agreement vifith the equal-mass, non- 
spinning numerical data, since the agreement is already 
at the level of the numerical error. Thus, in the following, 
we shall use the values of 05(1/4) and Wpoie(l/4) based 
on Fig. 2, obtained by setting to zero all the other EOB- 
dynamics adjustable parameters in Table I. We will leave 
a comprehensive study of the other EOB-dynamics ad- 
justable parameters to future work when highly accurate 
numerical merger waveforms of unequal-mass black-hole 
binaries become available. 

We shall now discuss the EOB model with reference 
values 05(1/4) = 6.344 and Wpoie(l/4) = 0.85, and 
tune the EOB-waveform adjustable parameters. We 
shall comment at the end of this section on the results 
when the other reference values 05(1/4) = 4.19 and 
i'poio(l/4) — > 00 are used. In Fig. 5 we compare the 
numerical and EOB /122 amplitudes with and without 
including NQC terms. The agreement of the numerical 
amplitude with the EOB amplitude of Eq. (16a) without 
NQC terms, which uses the resummation procedure of 
Ref. [17], is rather remarkable. The relative difference at 
the peak is only ~ 1.5%, and the EOB peak amphtude 
occurs only ^ 6M before the numerical peak amplitude. 



We notice that this excellent agreement is due to the 
presence in P22 of test-particle corrections through 5PN 
order. Were the test-particle corrections through 4PN or 
5PN orders not included, the disagreement at the peak 
would become 4.9% and 11.3%, respectively.^ 

Figure 5 also shows the EOB amplitudes of Eq. (16a) 
when the Pade-resummations P4 and P| of P22 suggested 
in Ref. [17] are applied. In these cases, the EOB peak am- 
plitude almost coincides in time with the numerical peak 
amplitude, but the relative difference in the value of the 
peak amplitude is rather large. However, those large dif- 
ferences may be resolved if the resummed version of the 
GW energy flux [17] consistent with the resummed him 
were used. Figure 5 also contains the EOB ampli- 
tude with NQC terms as suggested in Refs. [30, 31] [see 
Eq. (12) in Ref. [31]]. The relative difference with the 
numerical amplitude is ~ 20% at the peak. It is rather 
interesting to observe, as pointed out in Ref. [31], that 
by aligning the numerical and EOB waveforms at low 
frequency, we find that the peak of the numerical h22 co- 
incides with the peak of the EOB orbital frequency. Here, 
to improve the amplitude agreement during plunge and 
merger, we include the NQC corrections of Eq. (18). We 
fix two of the adjustable parameters, oj^^ and 03^^, by 
requiring that a local extremum of the EOB /i22 ampli- 
tude occurs at the same time as the peak of the EOB 
orbital frequency (i.e., the EOB light-ring), and that the 
EOB amplitude at the peak coincides with the numerical 
amplitude at the peak. In fact, we expect that in the near 
future, the peak of the numerical h22 will be able to be 
predicted by numerical relativity with high accuracy for 
several mass ratios. Thus, the peak can be fit with a poly- 
nomial in ly. (Preliminary studies which use results from 
Ref. [28] confirm this expectation.) The other two ad- 
justable parameters, i.e., 03^^ and 04^^, are calibrated to 
the numerical results to further reduce the disagreement. 
Specifically, we do a two-parameter least-square-fit of the 
ratio of the numerical RWZ and EOB /i£,„ on Eq. (18) in 
which o^^^ and al^^ are fixed as functions of 03^^ and 04^ 
by the requirements described above. We notice that the 
strategy of improving the amplitude agreement followed 
in this paper might change in the future, when accu- 
rate numerical unequal-mass black-hole binary inspiral- 
merger-ringdown waveforms become available. A smaller 
number of adjustable parameters might suffice if more re- 
quirements on the EOB model itself can be imposed or if 
a different matching procedure, such as the one suggested 
in Ref. [51], is employed. 



^ In Ref. [17] (sec Fig. 10 therein and discussion around it) the 
authors pointed out that the difference between /122 amplitudes 
computed with the test-particle corrections through 3PN, 4PN or 
5PN orders, differ only by a few percent. However, this statement 
was obtained for circular orbits until the LSO frequency Mf! = 
0.097. Our Fig. 5 extends beyond that frequency (the latter 
corresponds to t = 3914A/ in the figure). 
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FIG. 6: Comparison of numerical waveform to EOB waveform 
with a5(l/4) = 6.344 and Upoic(l/4) = 0.85, i.e. the same 
model used in Fig. 5. The top panels show the real part of 
numerical and EOB /i22, the bottom panels show amplitude 
and phase differences between them. The left panels show 
times t — to 3900M, and the right panels show times t — 
3900 to t = 4070M on a different vertical scale. 



Comparing the gravitational-virave modes hiri 
equal-mass coalescing black-hole binaries 



of 



In this section, we focus on the model whose EOB- 
dynamics and EOB-waveform adjustable parameters 
were calibrated to numerical RWZ h22 and NP ^f^^ in 
Sec. IIIB. Using this EOB model, we generate the GW 
modes /i22, ^32 and /144 and compare them to the corre- 
sponding numerical modes. We choose these three modes 
because they are the most dominant ones for an equal- 
mass, non-spinning black- hole binary. 

In Fig. 6, we show the numerical and EOB mode /i22 
aligned with the procedure of Sec. Ill A. Using the ref- 
erence values 05(1/4) = 6.344 and t'poic(l/4) = 0.85, 
we find that the best phase and amplitude agreement 
is obtained when the matching occurs at an interval of 
At^Lch = 3.0M ended at t^Ztch = 3942.5M, i.e., at the 
peak of Mn, with a^'(l/4) = -2.23 and a2''(l/4) = 
31.93, a3"'(l/4) = 3.66 and a4""(l/4) = -10.85. The 
phase difference is strictly within ±0.02 radians until 
the merger, i.e., the peak of /122, which happens at 
t ~ 3942. 5M (early numerical data contaminated by junk 
radiation was discarded until t = 200M). The relative 
amplitude difference is also within ±0.02 in this range. 
The phase difference becomes 0.04 radians a.t t — 3962M; 
before a rather large error starts contaminating the nu- 
merical /122 • A more careful tuning on the EOB-waveform 
adjustable parameters could further improve the phase 
agreement. However, we do not think it is worthwhile 
to improve the agreement at this point since we are only 
examining the equal-mass case. Note that the relative 
amplitude difference becomes ^ 7% a.t t — 3962 Af, and 
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FIG. 7: Comparison between EOB /122 and the numerical 
"tl^ . The top four panels show the real part of the waveform, 
on a linear and logarithmic j/-axis. The bottom two panels 
show the phase difference (in radians) and the fractional am- 
plitude difference between the two waveforms. The left panels 
show times t = to 3900M, and the right panels show times 
t = 3900 to t = 4070Af with different vertical scales. (The 
quantities in the lower left panel have been smoothed; the 
grey data in the background of that panel presents the raw 
data.) This figure uses the same EOB model as Figs. 5 and 6, 
namely a5(l/4) = 6.344 and 7;poie(l/4) = 0.85. 



increases during the ringdown. 

The numerical GW strain /122 plotted in Fig. 6 is com- 
puted using RWZ wave extraction. During the ringdown, 
this waveform is noisier than the extracted NP scalar '^^^ 
(see the Appendix). Therefore, in Fig. 7, we compare the 
numerically extracted '^^^ with the second time deriva- 
tive /i22 of the EOB waveform. Overall, the agreement 
is much better than for the comparison of /122 in Fig. 6. 
Phase and relative amplitude differences are smaller than 
0.002 during most of the inspiral, and remain smaller 
than 0.01 up to f = 3920M. In the interval around 
merger, t = 3930M to 3960A/, the agreement is slightly 
worse than in Fig. 6; the disagreement in this region is 
caused by the differences between the inspiral EOB /122 
and numerical NP 'I'l^ frequencies, as discussed at the 
beginning of Sec. III. 

In the ringdown region, t > 3960 A/, Fig. 7 shows ex- 
cellent agreement, and this agreement persists until late 
times. In contrast to the h comparison shown in Fig. 6, 
in Fig. 7 both phase and amplitude differences remain 
bounded; during the ringdown, the phase difference be- 



tween EOB /122 and '^^■^ oscillates around 0.08 radians. 
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FIG. 8: We show the amphtude and frequency of the numer- 
ical and BOB mode /i22, the EOB orbital frequency and the 
frequency of the numerical mode ^'4^ . The vertical line marks 
the peak of the EOB amplitude and orbital frequency. 



and the amplitude differs by about 8%. Apart from 
small oscillations likely caused by gauge effects (see the 
Appendix), A(/) remains constant to an excellent degree 
during about 9 ringdown oscillations, i.e. during an accu- 
mulated phase of about 56 radians. If the quasi- normal 
mode frequency used in the EOB ringdown waveform 
were different from the numerical ringdown frequency by 
as little as 0.1%, a linearly accumulating phase-difference 
of ~ 0.056 radians would accumulate, which would be 
clearly noticeable in the lower right panel of Fig. 7. Thus, 
we find agreement at the 0.1% level between the nu- 
merical quasi-normal mode frequency and the prediction 
based on final mass and spin of the numerical simulation. 

In Fig. 8, we compare the amplitude and frequency 
of numerical and EOB /i22 waveforms together with the 
orbital frequency of the EOB model. The peak of the 
latter is close to the EOB light ring, and is aligned with 
both the EOB and numerical /122 amplitudes (as required 
by our choice of a'j*^^ and 0-2^^) ■ During the ringdown, 
the frequency computed from the numerical ft.22 shows 
increasingly large oscillations. We also plot the frequency 
computed from the numerical mode. This frequency 
shows much smaller, and bounded, oscillations deep into 
the ringdown regime. 

Having constructed our EOB waveform purely by con- 
sidering the (2,2) mode, we now discuss agreement be- 
tween higher modes of the EOB model and the numerical 
simulation. Figure 9 shows phase and amplitude differ- 
ences for the two next largest modes, the (4,4) and the 
(3,2) mode. The EOB model is identical to the one that 
has been calibrated to agree with the (2,2) mode, and the 
parameters a^*^^ and aj*", which appear in Eq. (18) to 
correct the amplitude of the higher-order modes for non- 




FIG. 9: Upper panel: Amplitude and phase differences of nu- 
merical and EOB mode /132 over the inspiral range. Lower 
panel: Amplitude and phase differences of numerical and 
EOB mode /i44 over the inspiral range. 



quasi-circular motion, are set to zero. But although the 
EOB model has not been calibrated in any way to match 
the higher-order modes, the agreement between numeri- 
cal and EOB waveforms shown in Figure 9 is rather good 
for t < 3700M. In fact, the differences between EOB and 
NR modes are comparable to the estimated numerical er- 
rors in these modes (as estimated by convergence tests 
between different numerical resolutions, and the compar- 
ison between the numerical h and 4*4 waveforms which 
are presented in the Appendix. Around t w 3700M, the 
numerical (3,2) and (4,4) modes begin to show additional 
features, which we believe are unphysical, and are de- 
scribed in more detail in the Appendix. These features 
prevent a meaningful comparison of the (3,2) and (4,4) 
modes at later times. 

Figure 10 shows amplitude and frequency of the (4,4) 
and (3,2) modes for both the EOB model and the nu- 
merical simulation for the last few hundred M of inspiral. 
This figure begins approximately where the NR-EOB dif- 
ferences in Fig. 9 exceed the vertical scale of that figure. 
The EOB amplitude and phase follow roughly the aver- 
age of numerical results, which show oscillations resulting 
from numerical errors. At earlier times, the EOB and NR 
amplitudes, phase and frequencies track each other very 
closely, as can be seen from Fig. 9. Please compare also 
with Fig. 8 which plots the frequencies for the (2,2) mode. 

Finally, in Fig. 11, we show the numerical and EOB 
mode /i22 using the reference values 05(1/4) = 4.19 and 
"^poic — *■ 00. In this case we find that the best phase 
and amplitude agreement is obtained when the matching 
occurs over a range of Ai^^^j^i^ = 2.2Af ended at the peak 
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FIG. 10: Upper panel: Amplitude of the numerical and EOB 
modes /i32 and /144. Lower panel: Frequency of the numerical 
and EOB modes /132 and /144. The EOB orbital frequency 
2MQ (4A/f2) is indistinguishable from the frequency of the 
h32 (/144) mode on the scale of this plot. 
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FIG. 11: Comparison of the numerical data to an EOB model 
with Upoic = 00. This figure is analogous to Fig. 6, but uses an 
EOB-model that was calibrated with the restriction iipoic = 00 
(parameters are given in the main text). Even without Wpoic, 
the inspiral can be matched equally well as in Fig. 6; during 
the ringdown, the phase differences are somewhat larger, but 
it is possible that refined tuning will reduce them further. 



of MQ, with a'l'-^{l/4) ^ -2.50 and a2''(l/4) = 35.43, 
4=^" (1/4) = 4.91 and 04=^=^(1/4) = -32.40. Comparing 
the result with that of Fig. 6, we notice that the phase 
and amplitude differences are only slightly worse than 
the reference model of Fig. 6, but still within numerical 
error. 



D. Impact on data analysis 

Using the EOB model with reference values 05(1/4) = 
6.344 and ?;poie(l/4) = 0.85, we now quantify the dis- 
agreement between numerical and EOB waveforms by 
calculating their maximized overlaps which are impor- 
tant for analysis of data [52] from GW detectors. Here 
we restrict ourselves to the dominant mode /i22- Given 
two time-domain waveforms hi{t) and /i2(^; ^Oj </'o) gen- 
erated with the same binary parameters, the maximized 
overlap, otherwise known as a fitting factor (FF), is given 
explicitly by [27] 

FF ^ max (hi^h^jto^M) (34) 

to,0o ^ {hi, hi) (/l2(io, 0o), /i2(to, 0o)) 

where 

(/^i,M^4Re^°°^i^MZ)#. (25) 

Here hi{f) is the Fourier transform of hi{t), and Sh{f) 
is the detector's power spectral density. We compute the 



FFs for binary black holes with total mass 30-150A'/q, 
using LIGO, Enhanced LIGO and Advanced LIGO noise 
curves ^, and find in all cases FFs larger than 0.999. 
Note that the FFs are computed maximizing over time 
of arrival and initial phase, but not over the binary pa- 
rameters. We note that FF > 0.999 gives a mismatch 
e = 1 — FF between the numerical and the analytical 
/i22 of cnr-eob < 0.001. For the noise curves of LIGO, 
Enhanced LIGO and Advanced LIGO, we find that the 
mismatch between all extrapolated numerical waveforms 
h is less than 0.0001 for black-hole binaries with a total 
mass of 3O-I5OM0. If we take this mismatch as an es- 
timate of the difference between the numerical and the 
exact physical waveforms, we have Eo-nr < 0.0001. The 
mismatch between the exact and the analytical /i22 is 
therefore Ce-EOB < 0.0017. This mismatch is smaller 
than the bound 0.005 presented in Ref. [8], and therefore 
our EOB model is sufficiently accurate for GW detection 
in LIGO, Enhanced LIGO, and Advanced LIGO. 



® For LIGO, we use the analytic fit to the LIGO design power 
spectral density given in Ref. [53]; for Enhanced LIGO, we 
use the power spectral density given at http://www.ligo. 
caltech.edu/~raiia/NoiseData/S6/DCnois6.txt; for Advanced 
LIGO, we use the broadband configuration power spectral den- 
sity given at http://www.ligo.caltech.edu/advLIGO/scripts/ 
ref _des . shtml 
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E. Unequal mass inspiraling binary black holes 



As a check of the robustness of our EOB model cal- 
ibrated to numerical waveforms of equal-mass black- 
hole binaries, we extend the model to a set of unequal- 
mass black-hole binaries by comparing numerical and 
EOB inspiraling waveforms for mass ratios 2:1 and 
3:1. These simulations were performed with the Caltech- 
Cornell SpEC code, last about eight orbits and have 
phase errors similar to the equal mass simulation dis- 
cussed so far. Details of simulations will be published 
separately. We explore here the possibility of setting 
0-5 (v) ~ v\q with Ao constant and let Wpoio depend on the 
mass ratio. Indeed, in the test-particle limit we expect'' 
from l/v^. Wpoio(O) = l/-\/3 = 0.57735, whereas in the 
equal-mass case we find Wpoic(l/4) = 0.85. As a prelim- 
inary study, we do not perform a comprehensive search 
over the Ao-Wpoio parameter space for unequal-mass bina- 
ries, as we did for equal-mass binaries in Sec. IIIB. We 
fix Ao to our reference value 25.375 and tune v-po\c{v) to 
require phase differences on the order of the numerical 
error. 



In Figs. 12 and 13, we compare the numerical and EOB 
waveforms and their amplitude and phase differences 
for binaries with mass ratios q ~ mi:m2 of 2:1 and 3:1. 
The alignment procedure of Sec. Ill A was used with 
ti = 310M and ^2 = 930M. The figures also show the nu- 
merical phase error obtained from runs with two different 
numerical resolutions. The specifics of these numerical 
runs will be published separately. In the case of mass ra- 
tios q = 2:1 and 3:1, we find that by tuning Wpoie(i^), the 
difference between numerical and EOB waveforms can 
be reduced to values smaller than the numerical error. 
The best values of v^oie we find are Vpoic = 0.76 ± 0.01 
for mass ratio 2:1, and Wpoie = 0.70 ± 0.01 for mass 
ratio 3:1. Choosing parameters outside this range re- 
sults in differences between numerical and EOB wave- 
forms that are at least twice the numerical error. Com- 
bining Wpoio values for mass ratios 1:1, 1:2, 1:3, and the 
test-particle limit, we find a least-square fitting formula 
v^o\c{v) = 0.57 - 0.65(±0.35):/ + 7.Q{±l.b)v'^. 

Finally, we observe that the phase and amplitude dif- 
ferences between numerical and EOB waveforms can be 
reduced to values smaller than the numerical error, if we 
choose the EOB reference model of Sec. Ill B where we 
set Vpo\e oo and let az(v) = v{\q + Xiv)- In par- 
ticular, calibrating the mass ratio 2:1 and 3:1, we find 
a^{v) = i/[-7.3(±0.1) + 95.6(±0.3)i^]. These EOB mod- 
els agree with the numerical data as well as the EOB 
models shown in Figs. 12 and 13. 




^ We note that with our choice of the GW energy flux (factor- 
ized logarithms and vi^o = 1), the best fit to numerical flux has 
''polcCO) = 0.57, which differs slightly 
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FIG. 12: EOB-NR comparison for a BH binary with 
mass ratio 2:1. The upper panel shows the numerical 
and EOB mode , and the lower panel shows phase and 
amplitude differences between EOB and numerical run. The 
dashed brown line is the estimated phase-error of the numeri- 
cal simulation, obtained as the difference between simulations 
at high resolution 'N6' and lower resolution 'N5'. 



IV. CONCLUSIONS 

In this paper, building upon recent, successful re- 
sults [17, 22, 26-31] of the EOB formalism [9, 10, 13- 
15], we have concentrated on the EOB model denoted 
""^Fi/Hi in Ref. [22], with adjustable EOB-dynamics 
and EOB-waveform parameters defined in Table I. We 
have calibrated this EOB model to a very accurate nu- 
merical simulation of an equal-mass non-spinning binary 
black-hole coalescence [23] . 

When comparing EOB and numerical waveforms, or 
when comparing numerical waveforms with each other, 
we determine the arbitrary time offset and phase off- 
set between the waveforms by minimizing the phase dif- 
ferences between the waveforms over a time interval of 
^ lOOOM at low frequency, where the PN-based EOB 
waveforms are expected to be most accurate [22]. Com- 
pared to aligning waveforms at a particular time or fre- 
quency, this procedure is less sensitive to numerical noise 
and residual eccentricity. 

Among the EOB-dynamics adjustable parameters 
{a5(l/4),Wpoie(l/4),aJ|(l/4),a;^J^(l/4),^8}, the pa- 
rameters 05(1/4) and Wpoic(l/4) have the largest effect 
upon the long inspiral phase. Thus, we set {app(l/4) = 
0, a;^j^(l/4) = 0, ^8 = 0} in our EOB model, and we con- 
sidered the phase difference between the numerical and 
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FIG. 13: EOB-NR comparison for a BH binary with 
mass ratio 3 : 1. The upper panel shows the numerical 
and EOB mode ^4^, and the lower panel shows phase and 
amplitude differences between EOB and numerical run. The 
dashed brown line is the estimated phase-error of the numeri- 
cal simulation, obtained as the difference between simulations 
at high resolution 'N6' and lower resolution 'N5'. 



EOB as a function of the parameters 05(1/4) and 
Wpoic(l/4)- This phase difference increases with time, so 
we have sought parameters for which this phase difference 
remains small for as long a time as possible. We found re- 
gions of the (05(1/4), Wpoio( 1/4)) parameter space where 
this phase difference is less than 0.02 radians either until 
t = 282M or until t = 42M (red curves) before the time 
when the numerical /122 reaches its maximum amplitude 
(see blue and red curves in Fig. 2), respectively. 

Moreover, building on Refs. [10, 50], we have found 
that the EOB-dynamics adjustable parameters entering 
the GW energy flux cannot modify the phase of the EOB 
^'4^ during the plunge and close to merger. This is be- 
cause any modification of the GW energy flux beyond 
the LSO has negligible effect on the phasing, as the evo- 
lution is driven mostly by the conservative part of the 
dynamics. We also found that As is strongly degener- 
ate with 05(1/4) until almost lOOA/ before merger, and 
that app(l/4) and o^^(l/4) have an almost equivalent 
effect on the phasing, except for a minus sign and a dif- 
ferent scaling. Overall, for the equal-mass non-spinning 
case, we have found that the EOB-adjustable parameters 



{oj^^(l/4), a^^{l/A), Ag,} have a minor effect in reducing 
the phase and amplitude differences between the EOB 
model and the numerical simulation (see also Fig. 4). To 
achieve differences on the order of the numerical error, 
we can restrict ourselves to the EOB parameter space 



with {a^|(l/4) = 0,0^^,(1/4) = 0, ^8 = 0}. 

Furthermore, using our alignment procedure, we have 
found that the peak of the numerical ^,22 coincides with 
the peak of the EOB orbital frequency, confirming what 
was pointed out in Ref. [31]. As in Rcf. [31], wc re- 
quire that the EOB dominant mode /122 peaks at the 
maximum of the EOB orbital frequency (i.e., the EOB 
light-ring). We also require that the EOB amplitude at 
the peak coincides with the numerical amplitude at the 
peak. In fact, we expect that in the near future, the 
peak of the numerical /122 will be able to be predicted by 
numerical relativity with high accuracy for several mass 
ratios. Thus, the peak can be fit with a polynomial in v. 
(Preliminary studies which use results from Ref. [28] con- 
firm this expectation.) These requirements determine the 
EOB-waveform parameters a'^^^l/A) and a2''(l/4). To 
further improved the agreement close to merger, we then 
tune a^^^(l/4), 04^^(1/4), and At^^^^^{l / A) , so that the 
phase and amplitude differences between the EOB and 
numerical /122 ai'c minimized. In particular, we found 
that this happens if ^^^^^^^^1^(1/4) is chosen to be around 
3A/ (while imatch(l/4) is fixed at the maximum of the 
EOB orbital frequency Mfi). For the EOB reference 
model with 05(1/4) = 6.344 and Upoic = 0.85, we have 
found that the phase and amplitude differences between 
EOB and numerical /122 waveforms are 0.02 radians and 
2%, respectively, until 20M before merger, and 0.04 radi- 
ans and 7%, respectively, during merger and early ring- 
down, until the numerical /i22 starts to be affected by nu- 
merical oscillations (see Fig. 6). These agreements were 
obtained by comparing EOB and numerical values of /i22 , 
the latter having been extracted from the RWZ scalars. 
We also compared the EOB and numerical 'I'l^. In this 
case, the agreement is even better during the long inspiral 
and through the late ringdown, with phase and amplitude 
disagreements of 0.02 radians and 2% until 20A/ before 
merger, and 0.08 radians and 8%, respectively, during 
merger and ringdown (see Fig. 7). However, around the 
transition between plunge and ringdown, the EOB /i22 
has some oscillations because the EOB resummation pro- 
vides us with /i22, whereas when taking time derivatives 
of h22 non-resummed higher order PN terms are gener- 
ated, spoiling in part the agreement of /i22. 

Quite interestingly, we have found that phase and am- 
plitude differences between EOB and numerical wave- 
forms can also be reduced to numerical errors, at least 
during the inspiral, if we let Wpoic — *■ 00 and calibrate the 
coefficients Aq and Ai in 05(1^) = i/(Ao + Aii/), see Fig. 11. 

For data analysis purposes, wc have also computed 
the maximized overlaps or fitting factors (EEs) between 
the EOB reference model with 05(1/4) = 6.344 and 
fpoic ^ 0.85 and numerical /i22. We maximized only over 
the initial phase and time of arrival. We have found that 
for black-hole binaries with total mass 3O-15OA'/0, us- 
ing LIGO, Enhanced LIGO and Advanced LIGO noise 
curves, the EEs are larger than 0.999. We have also com- 
puted the EEs between values of numerical h22 that were 
computed in slightly different ways (e.g. different nu- 
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merical resolutions, different extraction procedures), and 
have estimated tlie misniatcli between tlie exact and EOB 
/i22- We have concluded, in the spirit of Ref. [8], that our 
analytical /i22 satisfies the requirements of detection with 
LIGO, Enhanced LIGO and Advanced LIGO. 

Finally, to test the robustness of the EOB model, we 
have also compared it to a few equal-mass subdominant 
modes {£,m), notably (4,4) and (3,2), and to the domi- 
nant mode (2, 2) of a set of unequal-mass inspiraling bi- 
naries. Without changing the EOB-dynamics adjustable 
parameters, we have found that, in the equal-mass case, 
the phase and amplitude differences of EOB and numer- 
ical /144 and /132 are within the numerical errors through- 
out the inspiral (see Figs. 9 and 10). Furthermore, in the 
unequal-mass case, we have found that we can reduce the 
phase difference of the EOB and numerical '^^^ of inspi- 
raling binaries of mass ratios 2:1 and 3:1 on the order of 
the numerical error (see Figs. 12 and 13). This can be 
obtained either (i) by setting 05 (i^) = vXq with Aq fixed 
by the equal-mass case, and calibrating fpoie(j^), or (ii) 
by letting Wpoic ^ 00 and calibrating 05 (i^). 

In the near future, we plan to compare the non- 
spinning EOB model defined in this paper to a larger 
set of accurate numerical simulations of black-hole bi- 
nary coalescences (for both equal and unequal-mass 
binaries), and complete the tuning of all the EOB- 
dynamics and -waveform adjustable parameters. In par- 
ticular, we expect to improve the EOB plunge-mergcr- 
ringdown matching either by reducing the number of 
EOB-waveform adjustable parameters or by employing 
different matching procedures or GW energy fluxes. 

While polishing this manuscript for publication, an in- 
dependent calibration of the EOB model which uses the 
equal-mass binary black-hole data of the Caltech-Gornell 
collaboration employed in this paper and made public on 
January 20, 2009 appeared on the archives [54]. 



Acknowledgments 

We thank Oliver Rinne for his work on implement- 
ing Regge-Wlieeler-Zerilli wave extraction, and Fan Zang 
for extrapolating waveforms to infinite extraction radius. 
We also thank Emanuele Berti and Evan Ochsner for 
useful discussions, and Emanuele Berti for providing us 
with the quasi-normal mode frequencies and decay times 
used in this paper. A.B. and Y.P. acknowledge support 
from NSF Grant No. PHY-0603762. L.B., L.K., H.P., 
and M.S. are supported in part by grants from the Sher- 
man Fairchild Foundation to Caltech and Cornell, and 
from the Brinson Foundation to Caltech; by NSF Grants 
No. PHY-0601459, No. PHY-0652995, and No. DMS- 
0553302 at Caltech; by NSF Grants No. PHY-0652952, 
No. DMS-0553677, and No. PHY-0652929 at Cornell. 



APPENDIX: COMPARING DIFFERENT 
METHODS OF COMPUTING hirrv 

The analysis in Sec. Ill relies to some extent on the GW 
strain h extracted from the numerical simulation. Earlier 
papers describing generation of the numerical data [21- 
23] focused on the behavior of the NP scalar ^'4, and 
performed comparisons to PN theory based on the nu- 
merical 5*4. 

We have two means of computing a GW strain h from 
the numerical simulations. The first is a double time 
integration of 5*4, exploiting the relation 

*4 = h. (A.l) 

[Note that throughout this Appendix, we suppress in- 
dices £jn denoting the components of the decomposition 
into spin- weighted spherical harmonics. Thus, Eq. (A.l) 
is meant to apply to each complex component {£,m)]. 
For each time integration [and each mode {i,m)], a 
complex integration constant needs to be determined. 
These constants are fixed with the procedure described 
in Sec. II of Ref. [22], in which a certain functional of 
temporal variations of the amplitude of the integrated 
data is minimized. The minimization is performed over 
25 separate integration intervals [ti , t2] with ti /AI ~ 
1000, 1100,..., 1400 and ta/M = 2600, 2700, 3000. 
We then compute the time average of these 25 integrated 
waveforms, and we use this time average, which we de- 
note as // 4'4, as the GW strain. Note that we perform 
the above operations on the numerical ^'4 data after it 
has been extrapolated to infinite extraction radius. 

Our second means of extracting a GW strain is us- 
ing the RWZ equations [44, 45] generalized to arbitrary 
spherically symmetric coordinates, as formulated by Sar- 
bach & Tiglio [46] . An advantage to the Sarbach & Tiglio 
formalism in contrast to the more widely-used Zerilli- 
Moncrief formalism ( [55] and references therein) is that 
in the former case, the GW strain is obtained directly 
from the gauge-invariant RWZ scalars (at leading order in 
the inverse radius), without any time integration. With 
Oliver Rinne, we have implemented the Sarbach & Tiglio 
formalism for a Minkowsi background in standard coor- 
dinates in the Caltech-Gornell spectral code [47]. From 
the RWZ scalars (extracted at finite radii), we compute 
the GW strain and then extrapolate to infinite extraction 
radius in order to obtain the final waveform /irwz- 

In order to gain insight into the accuracy and relia- 
bility of the computed GW strain, we explore the differ- 
ences between waveforms extracted with either technique 
(see also [56] for a similar comparison). Figure 14 shows 
the real part of the numerical (2, 2) mode. On the scale 
of the full waveform, no disagreement between /irwz and 
// ^'4 is visible. However, the lower two panels of Fig. 14 
show differences between /irwz and 'I'4 deep in the ring- 
down phase: While ^'4 continues to decay exponentially 
through many orders of magnitude, /irwz exhibits notice- 
able deviation from a pure exponential decay at about a 
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FIG. 14: The (/, m) = (2, 2) mode of the numerical waveform. 



tenth of peak amplitude. Decay of /irwz stops completely 
at about one per cent of peak amplitude. 

We suspect that this unexpected behavior is caused by 
gauge effects: All simulations in the numerical relativity 
community are performed using gauges in which the co- 
ordinates dynamically respond to the changing geometry, 
so as to avoid pathologies such as coordinate singulari- 
ties. Ideally, the procedures used to extract gravitational 
radiation from the simulations should be gauge invari- 
ant, so that the choice of gauge used in the simulation 
is irrelevant. In practice, however, wave extraction tech- 
niques are not perfect. For example, the RWZ technique 
is gauge invariant only to first order in perturbation the- 
ory about fixed background coordinates. Likewise, the 
NP technique is strictly gauge-invariant only if applied 
at future null infinity, rather than at a finite distance 
from the source. Gauge effects are expected to manifest 
themselves differently in NP and RWZ wave extraction 
techniques, so by comparing the results of these two ex- 
traction techniques, we can get a handle on the size of 
our uncertainties that arise from gauge effects. 

Therefore, we will examine the differences between the 
numerical /irwz and 5*4. Using (A.l), we can compute 
a meaningful difference in two ways. The first way is to 
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FIG. 15: Phase and relative amplitude difference between the 
{I, m) — (2, 2) modes of the RWZ waveform /irwz and NP 
scalar ^4, [see Eqs. (A.2)-(A.5)]. The right panel shows an 
enlargement of merger and ringdown, with the dotted vertical 
lines indicating time of maximum of |^'4|, and where |vl'4| has 
decayed to 10% and 1% of the maximal value. (The blue lines 
are smoothed; the grey data in the background represents the 
unsmoothed data.) 



differentiate /irwz twice and compute 

A0NP = arg(*4) - arg(/iRwz) 

A^NP ^ |^4| - l/tRWzl 
^ ~ (|W| + |*4|)/2' 



(A.2) 
(A.3) 



The subscript 'NP' indicates that the comparison is made 
on the level of the NP scalars, i.e. ^"4 appears undiffer- 
entiated on the right-hand-sides. The second way is to 
time-integrate 5*4 and to calculate 



A0RWZ arg(// ^"4) - arg(/iRwz) 
A^RWZ _ I// ^4| - I^RWzl 
A ~ (|W| + |//^4|)/2- 



(A.4) 
(A.5) 



The results of these comparisons are presented in Fig. 15. 
An examination of this figure reveals several properties 
of the extracted ^'4 and /irwz waveforms. First, we note 
that during the inspiral and merger (up to t < 3960M, 
that is 18M after the peak of /irwz), the RWZ and NP 
waveforms agree to better than 0.02 radians. A(/)np con- 
tains more noise because noise is amplified by the double 
time differentiation to compute /i, and because ^"4 is con- 
taminated by junk-radiation from the initial data up to 
time t w lOOOAf . The blue lines in this plot have been 
smoothed (by convolution with a Gaussian of width 5M) 
to reduce the effect of noise due to junk radiation. (The 
grey data in the background of Fig. 15 shows the un- 
smoothed A$Np)- In contrast, A(/)rwz does not show 
similar high frequency noise (the red dashed curves in 
Fig. 15 are not smoothed). Integration naturally smooths 
noise and apparently, the RWZ wave extraction is less 
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FIG. 16: The (/,m) = (4, 4) mode of the numerical waveform. 



susceptible to the noise introduced by junk radiation. 
Unfortunately, because of an imperfect choice of integra- 
tion constants for the time integration, //4'4 docs not 
precisely oscillate around zero at all times. This results 
in oscillations of A0rwz and AArwz/^ during the in- 
spiral; the frequency of these oscillations coincides with 
the GW frequency. The choice of integration constants, 
however, is good enough to confine these oscillations to 
less than about 0.02 radians in phase and 0.5 per cent in 
amplitude during the inspiral. 

Around merger, differences of the wave strain, i.e. 
and AArwz/^, begin to grow, and during ring- 
down this growth accelerates. This large disagreement 
is caused by two effects. The first effect is the contam- 
ination of /iRwz in the ringdown phase, presumably by 
gauge effects, as shown in Fig. 14. The second effect 
is related to the time integration used to obtain //^'4. 
During the inspiral phase, with an appropriate choice of 
integration constants the average value of // 4*4 is very 
nearly zero (see top left panel of Fig. 15). Thus, the in- 
spiral phase fixes all integration constants. When we now 
extend the integration through merger and ringdown, we 
find that / J ^"4 during ringdown has a contribution that 
grows linearly in time. Because the desired oscillatory 
part of // ^'4 decays exponentially, this linearly growing 
contribution contaminates arg J J ^1^4 to an increasing de- 
gree as time increases. The linearly growing contribution 
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FIG. 17; Phase and relative amplitude diflference between the 
(i, m) — (4, 4) modes of the RWZ waveform /lawz and NP 
scalar ^4, cf. Eqs. (A.2)-(A.5) The right panels shows an 
enlargement of merger and ringdown, with the dotted vertical 
line indicating the position of the maximum of |*4|. 



to // 5*4 is just barely visible in the top panel of Fig. 14; 
for the (3,2) and (4,4) modes discussed below, it will be 
much more obvious. 

When matching an analytical model waveform to nu- 
merical results, one must choose whether to match to 
^'4, // \l/4, or /iRwz, and we have just seen that these 
three numerical waveforms differ by systematic effects 
that arise from properties of the numerical simulation. 
Given Figs. 14 and 15, it appears that ^'4 is preferable 
over //^'4 because ^'4 lacks the low- frequency oscilla- 
tions during inspiral that are introduced in // 5*4 by time 
integration, and furthermore 5*4 lacks the linear drift dur- 
ing the ringdown. Similarly, ^'4 has an advantage over 
^Rwz because it has much cleaner behavior during ring- 
down (see Fig. 14). 

We now turn our attention to the next largest mode, 
{l,m) — (4,4), which is shown in Figure 16. Concen- 
trating on the top panel first, we see that J J ^4 agrees 
with /iRwz very well for a large fraction of the inspiral. 
However, for t < lOOOM and t > 3900M, // 1-4 contains 
contributions that grow linearly in time. Note that these 
contributions cannot be removed by a different choice of 
integration constants, because integration constants re- 
sult in addition of a Hnear term a + bt uniformly at all 
times. Hence, if the integration constants were changed 
to yield agreement for t < lOOOAf , the linearly growing 
discrepancy would appear at i > 1000 A/. The reason 
that the transition is around t ^ lOOOAf may be related 
to the so-called junk radiation that is present in numer- 
ical simulations, and arises because the initial data do 
not correspond precisely to a snapshot of an evolution. 
A small fraction of the outgoing junk radiation is re- 
flected when passing through the outer boundary. The 
reflected waves pass through the computational domain 
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(l,m) = (3,2) 
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FIG. 18: The (Z, m) — (3, 2) mode of the numerical waveform. 

at retarded time t « 100071/. While the reflected junk ra- 
diation is small, apparently it is sufficient to contaminate 
// ^'4, as seen in the top left panel of Fig. 16. 

Around merger, t sa 3950M, //^'4 picks up another 
linearly growing contribution which renders //^'4 basi- 
cally useless during merger and ringdown. This contam- 
ination might be related to oscillations in ^"4 and /irwz 
that become visible at t > 3750Af (see middle and lower 
panel of Fig. 16). It is presently unclear what causes these 
effects, but we conjecture that they are related to gauge 
effects that influence either our current wave-extraction 
procedure, or our current wave-extrapolation procedure. 
It is quite possible that a refined understanding of gauge 
effects will reduce these features in the future. 

Because of the apparent contamination of the wave- 



forms for early and late times, we will restrict the EOB- 
NR comparison of the higher order modes to the time in- 
terval 1000 < t/M < 3600. Figure 17 shows that within 
this interval, ^'4 and /irwz agree to better than 0.02 ra- 
dians in phase and 1% in amplitude. 

Finally, Figures 18 and 19 present an analogous com- 
parison for the (/, m) = (3, 2) mode. Qualitatively, these 
figures are similar to Figs. 16 and 17. Agreement be- 
tween ^4 and /iRwz is very good for the time interval 
1000 < t/M < 37GGAf, with the phases differing by less 
than 0.1 radians and the amplitudes by less than about 
1%. The larger disagreement might be due to the smaller 
RWZvs. - (I,m)=(3,2) 




1000 2000 3000 3600 3800 4000 



t/M t/M 

FIG. 19: Phase difference between the (l,m) = (3,2) modes 
of the RWZ waveform h and NP scalar *4, cf. Eqs. (A.2)- 
(A.5). The right panels shows an enlargement of merger and 
ringdown, with the dotted vertical line indicating the position 
of the maximum of |'I'4|. 



amplitude of the (3,2) mode of ^4 during the inspiral 
phase relative to the (4,4) mode. One potentially inter- 
esting difference between the (4,4) and (3,2) modes lies in 
the relative size of the variations in |ft.Rwz| and 1^*4 1 in the 
time range 3700 < t/M < 3900A/: For the (4,4) mode, 
variations in |^4| are clearly smaller than variations in 
I^RWzl (see Fig. 16. For the (3,2) mode this is reversed, 
with I ft. RWZ I showing somewhat smaller variations than 

1*41. 
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